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Tensor models generalize random matrix models in yielding a theory of dynamical triangulations 
in arbitrary dimensions. Colored tensor models have been shown to admit a 1/N expansion and 
^SJ ■ a continuum limit accessible analytically. In this paper we prove that these results extend to the 

most general tensor model for a single generic, i.e. non-symmetric, complex tensor. Colors appear 
, in this setting as a canonical book-keeping device and not as a fundamental feature. In the large A'^ 

■ limit, we exhibit a set of Virasoro constraints satisfied by the free energy and an infinite family of 
^ ■ multicritical behaviors with entropy exponents 7™ = 1 — 1/m. 
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I. INTRODUCTION 

JL . Matrix models are probability measures for random matrices AI of size N . In physics language, they come with 
a matrix action S{M). They can be divided in two broad categories. In the first category, that of invariant matrix 
(~| . models [l|, the full action has an expansion in terms of traces of powers of M (for Hermitian. or Tr(MM'l')" for 
general M) which ensures invariancc under U{N) transformations. The archetypes for this category are the ATr 
or ATr models, whose actions are 5 = TrM^ ATrM^ or S = TrM^ ATr M^. The perturbative expansion of 
^ ■ such models involves ribbon graphs dual to triangulated, or quadrangulated, Riemann surfaces. Hence (forgetting for 
I a brief moment the constructive issues) these models are statistical models of random discretized Riemann surfaces. In 
. the large N limit, planar surfaces dominate and furthermore undergo at some finite coupling a transition to continuous 
\0 ■ surfaces 0,11], known as the large volume, or continuum limit. Hence they provided until recently the only known 
, example of analytically controlled geometrogenesi^, i.e. the emergence of continuous geometries from discrete models, 

■ although restricted to two dimensions. Moreover invariant single or multi-matrix models can also probe the critical 
I behavior of two-dimensional statistical models on random geometries [ll-j^l- 

■ The second category of matrix models is that of matrix field theories, in which the interaction is invariant, but the 
quadratic part of the action is not. Since invariant Tr M" interactions are the matrix analogs of local interactions 
J ^"(x), matrix field theories are the analogs of ordinary quantum field theories, in which interactions are local but the 
propagator (inverse of the Laplacian or Dirac operator) is not. From this point of view invariant matrix models should 
be considered as ultralocal matrix field theories. Non-local propagators in field theory give birth to renormalization, 

■ hence to a flow of the couplings. Just as (t)\ is the archetype for ordinary renormalization, the archetype of matrix 
field theories is the Grosse-Wulkenhaar model in four dimensionsEI, or GW^ [l^l- This GW4 model improves on the 
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ordinary 04 model since it is asymptotically safe |14| , hence free of the old Landau ghost problem. 

Returning now to the important constructive question, let us recall that the constructive analysis of stable invariant 
matrix models is compatible with their 1/N expansion. Borel summability has been proved to hold uniformly in TV in 
the quartic case [l^. For higher degree stable interactions a straightforward generalization of the techniques of [l^ 
should also lead to uniform Borel-LeRoy summability of the appropriate order. The constructive analysis of matrix 



field theories is under way [17[ and expected to lead to a full construction of the C?W4 model in the near future. 

All these nice properties of matrix models stem from their \/N expansion [l8| . which states that planar graphs 
(dual to the sphere) govern their large limilE|. Planar graphs proliferate only exponentially in their number of 
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^ This term has appeared for the first time in [J which is however quite different from our approach. 

^ This critical behavior on random geometry is related to the one on fixed geometry through the KPZ correspondence [g| 4l2ll . 
^ The Grosse-Wulkenhaar model is a </<J^ model on the non-commutative Moyal space with a harmonic potential. It does not suffer from 
the UV/IR mixing and becomes a matrix field theory in the Moyal matrix base. 

* Through double scaling limits one can even to some extent treat the sum over sub-leading terms in the 1 /N expansion Il9l42l|l . 
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vertices and can be counted precisely through algebraic equations [22|, as they are related to trees |23l - l25| . This key 
feature underlies all the statistical mechanics applications of the invariant models. Renormalizability and asymptotic 
safety in the GW4 model also rely entirely on the dynamical analog of the expansion 

[l3,[Il[2i,[23- Indeed 

in such matrix field theories only planar graphs with a single external face look like matrix invariant terms at high 
energy, and they are also the only ones to require renormalization. 

Random matrices generalize in higher dimensions to random tensors psl - fsot (and [3ll - [33l | for more recent devel- 
opments), whose pcrturbativc expansion performs a sum over random higher dimensional triangulations. But until 
recently all nice aspects of matrix models listed above could not be generalized to tensors, as their 1/7V expansion 
was missing. The situation has changed with the discovery of colored (33-l36j rank D > 3 random tensor modcltQ. 
These models require D + 1 different pairs of conjugate tensors T', T*, equipped with a particular invariant canonical 
action of the type J2iT'f^ + AH.^o^' + ^IlZo'^'- Their perturbation theory supports a 1/N expansion ^MM, 
indexed by the degree, a positive integer which plays in higher dimensions the role of the genus but is not a topological 
invariant. Leading order graphs triangulate the Z?-dimensional sphere in any dimension (ssl. [39j. These graphs, bap- 
tized melonic [4]| . again proliferate only exponentially, as they map to colored {D + l)-ary trees (4ll.l4^. Like matrix 
models, these tensor models reach a continuum limit when the coup ling constant approaches its critical value. The 



corresponding entropy exponent is 7meions = 1/2 in any dimensions [41[. It is the analog of the string susceptibility 
exponent 7string = ^1/2 of the invariant matrix models for the universality class of pure 2d quantum gravity. 

Colored random tensors [4^ therefore gave the first theory of random geometries in three and more dimensions with 
analytically tractable geometrogenesis and the subject is rapidly expanding I44l450ll . Coupling of statistical mechanical 
systems to these random geometries in arbitrary dimension has been done in |48l . [sil . [52j , and results at all orders in 
1/N have been established for some restricted models [H^l (see also [s^l for some related developments). 

Obvious questions then arise. Do expansions also hold for uncolorcd models, i.e. with a single tensor? How 
can one build tensor models with interactions of arbitrary degree that still admit a 1/N expansion? What arc the 
tensor analogs of matrix field theories? 

A first important step towards answering the first two questions was taken in It was shown that integrating 
out all colored tensors but one in the initial colored model leads to an effective action for a single uncolored tensor 
which is a sum of effective invariant interactions whose internal structure can be unfolded in terms of colored graphs. 

In the present paper we return to these questions in greater detail. Like for matrices, we can distinguish invariant 
tensor models and tensor field theories. Invariant tensor models are those considered in this paper. They correspond to 
tensors with both quadratic part and interactions invariant under the external tensor product ®^U{N). We consider 
the most general invariant models for a random, complex tensor. It is important that this tensor is generic, that is 
without any symmetrization or antisymmetrization of its indices. Labeling these indices then provides exactly the 
same combinatorial tool that colors provide in the colored models. It allows us to 

• define their expansion, again organized according to the degree of the graphs, 

• prove it is dominated by melonic, colored graphs of spherical topology, 

• derive the continuum limit, whose entropy exponent is generically 7moions = 1/2, thus proving the universality 
of this continuum phas(|^, 

• extract a set of Virasoro constraints which hold in the large N limit, 

• find multicritical points, with entropy exponents 7m = 1 ~ (for m > 2 integers), which are the same as the 
ones of multicritical branched polymers |4^, [U, • This is the generalization of Q to tensors. 

These are the main results of this paper, and they are direct consequences of the universality of tensor invariant 
measures first derived in (56| . We stress that universality in this context only means that in the large N limit the 
tensors are distributed on a Gaussian. However, the Gaussian itself (i.e. its covariance) is not universal, but depends 
on the coupling constants of the model. Indeed, when the large N covariance becomes critical, the continuum limit 
is reached. Further, tuning the couplings appropriately, multicritical behaviors are observed, just like in invariant 
matrix models Q. 

We expect the constructive analysis of stable and symmetric invariant models not to pose any difficulty, as the 
necessary techniques have been in fact already developed for the quartic case [s^ in the slightly different context of 
group field theory [58{ . 



^ In D = 2 colors do not play the key role they play in three and more dimensions. This is because there is a natural composition rule on 
rank-2 tensors, namely matrix multiplication, and a single trace invariant at order n, namely TrM". In D > 3 there is no longer any 
multiplication law and there are many different invariants at order n. The colors become essential as a canonical device to keep track 
of their combinatories. Colored matrix models can of course still been defined and have been studied in [37j . 

^ It is analogous to the universality class of pure 2d quantum gravity which is obtained for most values of the coupling constants in 
one-matrix models [l|, Q. 
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Tensor field theories are the analogs of matrix field theories. They have tensor invariant interactions but a Laplacian- 
based propagator. Such a propagator again allows a renormalization group analysis. We do not consider this second 
category of models further in this paper, except to recall that uncolored renormalizable models of this type have 
been found for rank 3 and rank 4 tensors (sol . md^ . Again the renormalization in such models is entirely based on a 
dynamical version of the 1/iV expansion. One should explore their flows, phase transitions, critical exponents, gauge 
invariant extensions and constructive properties, as they seem a promising approach to the quantization of gravity in 
more than two dimensions j6l| . 

We follow the standard presentation of invariant matrix models in the large N limit, like in the well-known review 
to emphasize the new status of the field. In section HIl we define the generic models. In section Hill we consider 
their 1/N expansion, which is dominated by the melonic graphs and establish their continuum limit. In section ITVl 
we analyze the infinite family of multicritical points for these models. 



II. THE 1/iV EXPANSION OF INVARIANT TENSOR MODELS 
A. Tensor invariants and action 



The models we consider are based on complex tensors which have no symmetry between their indices. In order to 
write the most general action, one must first understand the invariants built from such tensors. It turns out that the 
analysis of these invariants automatically leads to a representation in colored graphs. 

Let Hi, . . . , Hjj be complex vector spaces of dimensions Ni, . . . , Njj. A rank D covariant tensor Tn-^^^ no can be 
seen as a collection of YiiLi complex numbers supplemented with the requirement of covariance under base change. 
We consider tensors T transforming under the external tensor product of fundamental representations of the unitary 
group <S)fLiU{Ni), that is each U{Ni) acts independently on its corresponding Hi. The complex conjugate tensor 
Tni...nD is then a rank D contravariant tensor. They transform as 

'^ai...an ~ ^ ^ C^aini ' ' ' V^dtid Tni...nn i ^ai...ar> ^ ^ ] Ua^no ' ' ' ^aini Tni...nD ■ (2-1) 

From now on wc will always denote the indices of the complex conjugated tensor with a bar. We will sometimes denote 
the D-uple of integers {ni, . . . ,nu) by ft and assume (unless otherwise specified) D > 3. Wc restrict to Hi = H, 
N, = N, for all i. 

Among the polynomial quantities one can build out of T and T we will deal in the sequel exclusively with trace in- 
variants. The trace invariants are built by contracting two by two covariant with contravariant indices in a polynomial 
in the tensor entries. We write trace invariants formally like 

Tr(r, T) = ^ n^»i^^i ^"1- ■ • ■ ' (2.2) 

where all indices are saturated. Note that a trace invariant has necessarily the same number of T and T. 

Trace invariants can be labeled by graphs with distinguished vertices. To draw the graph associated to a trace 
invariant we represent every T by a white vertex v and every T by a black vertex v. We promote the position of an 
index to a color: ni has color 1, 712 has color 2 and so on. The contraction of two indices n.^ and of tensors is 
represented by a line P = (v,v) connecting the corresponding two vertices. Lines inherit the color of the index, and 
always connect a black and a white vertex. Any trace invariant is then represented by a Z?-colorcd graph. 

Definition 1. A closed D-colored graph, or D-bubble, is a graph B = (V,£) with vertex set V and line set £ 
such that 

m V is bipartite, i.e. there exists a partition of the vertex set V = AU A, such that for any element I g £ , then 
I ~ {v, v} with V Cz A and v € A. Their cardinalities satisfy \V\ ~ 2|A| = 2|A|. 

• The line set is partitioned into D subsets E = E\ where is the subset of lines with color i, with |f | = |^| . 

• It is D-regular (all vertices are D-valent) with all lines incident to a given vertex having distinct colors. 

Some examples of trace invariants for rank 3 tensors are represented in figure [TJ The trace invariant associated to 
the graph B writes as 

D 

{n",n"}„,ijev v,v£B i=l l'' = {v,v)eB 
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where Z' runs over all the lines of color i of B. We call the (the product of delta functions encoding the index 

contractions of the observable associated to the graph B) the trace invariant operator with associated graph B [4^ . 
Trace invariant operators factor over the connected components of the graph. From now on we will always consider 
connected invariants, hence invariants associated to connected graphs in the above representation. We denote ^^^^ 
the set of D-colored, connected graphs with 2k distinguished vertices and r*^"^^ the set of all graphs with D colors. 

Of particular importance in the sequel are the subgraphs with two colors of a Z3-colored graph, called faces. We 
denote them F. For instance the graphs with 3 colors posses three type of faces, given by the subgraphs with lines 
of colors 12, 13 and 23. As every line belongs to exactly two faces (the lines of color 1 belong to a face 12 and a face 
13, etc.), the graphs with three colors can be represented as ribbon graphs. 

To every graph B with D colors we can associate a non- negative integer, its degree uj{B) [1^, H^, |4^. We recall 
its definition and properties in the appendix The main feature of the degree is that it provides a counting of the 
number of faces of a graph, thus for a graph with D colors and 2p vertices the total number of faces computes 

Taking into account that graphs with 3 colors are ribbon graphs, it is easy to see that in this case the degree reduces 
to the genus. In higher dimensions the degree provides a generalization of the genus. It is not a topological invariant, 
but it combines topological and combinatorial information about the graph. 

Going back to invariants one can build out of a complex tensor, we note that there exists a unique D-colored graph 
with two vertices, namely the graph in which all the lines connect the two vertices. We call it the D-dipole (denoted 
Bi) and its associated invariant is 

D 

Trs, (T, T) = ^ [n <5„,ft.] . (2.5) 

The most general invariant action for a non-symmetric tensor is therefore 

oo 

S{T, f) = ti Tre, (T, T) + ^ ^ te N~ro^^('^^ Tt:b{T, T) , (2.6) 

where (tjs) is the set of coupling constants associated to _D-bubbles and we singled out the quadratic part corresponding 
to Bi. In equation (|2.6p we have added a scaling in N for every trace invariant, proportional to its degree. As the 
degree is non-negative this scaling is a suppression of some invariants. We have included it because it simplifies some 
equations in the following, but we emphasize that this scaling is not required: as the reader can check all the results 
we present below generalize (albeit with some effort) also in its absence. Due to symmetry under relabeling of the 
black and white vertices, some couplings in (|2.6p are redundant. It is however more convenient to assign a distinct 
coupling constant to each graph with labeled vertices, and remember this redundancy only at the end. 

We will deal in this paper with the most general single-tensor model of rank D defined by the partition function 

Z{tB) = exp{-F{tt3)) = J dfdT exp (-N^-^ S{T,T)j . (2.7) 



B. Graph amplitudes 



The invariant observables are the trace invariants represented by /^-colored graphs. The Feynman graphs con- 
tributing to the expectation of an observable are obtained by Taylor expanding with respect to and evaluating 
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the Gaussian integral in terms of Wick contractions. A moment of reflection reveals that the Feynman graphs are 
made of effective vertices Tr b{T,T) (that is graphs B with colors 1, . . . ,D) connected by effective propagators (Wick 
contractions, pairings of T's and T"s). A Wick contraction of two tensor entries Ta^ . . ^^ and Tp^ . . ^^ with the quadratic 

part (|2.5[) consists in replacing them by js^D-if^ YiiLi ^atpt ■ The Wick contractions will be represented as dashed lines 
labeled by the fictitious color 0. Thus every dashed line of color in a Feynman graph identifies all the indices of the 
two vertices (one white corresponding to T and one black corresponding to T) it connects. An example of a Feynman 
graph is presented in figure [2l 

The Feynman graphs are therefore {D + l)-colored graphs Q. We reserve the notation B for the ZJ-colorcd graphs, 
and g for the {D + l)-colored graphs. A graph Q has two kinds of faces: those with colors i,j = 1, . . . , D, denoted 
Tij (which belong also to some D-bubble B) and those with colors 0, i, for i = 1, . . . , D, denoted J-Qi, which involve 
the lines of color in Q. 

The free energy has an expansion in closed, connected {D + l)-colorcd graphs. 

Fits)- E Hfir^^^^' ('-'^ 

where s(Q) is a symmetry factor and |p| is the number of effective vertices i.e. D-bubbles (subgraphs with colors 
1, . . . , D). We denote these _D-bubbles B^p^, with p = 1, . . . , \p\. The amplitude of a graph is 

^(^) = n*-<. E [n^"-"^"''<^''fen][ n r^n^«r,^r]- (2-9) 

An index is identified along the lines of color i in B(p') and along the dashed lines of color 0. We thus obtain a free 
sum per face of colors Oi, so that 

A{g) = ^^^(P) Jy{D-l)\p\-JshyJ:^u^iB^,^)-{D-l)\l°\ + J:^\J^„,\ _ j^Q-, 
J.l^'* I 

Noting that \ J^oi\ = \J- \ — J2p I-^(p)I' where denotes the total number of faces of G and |^(p)| the number of 
faces of the £)-bubble S(p), using (|2.4p for each K(p) and for Q (taking into account that Q has D + 1 colors) and 
noting that \l^\ = p, with p the half-number of vertices of Q, the amplitude of Q computes 

A{g) = iv^-iu^-(e) ^ (2.11) 

^1 

with ijj{Q) the degree of the graph Q. The 1/A^ expansion of the free energy writes 

F{ts) = N'' E (_^np^7V-(^"(^) . (2.12) 

The leading scaling with TV of the free energy is F{tis) ^ . In the rest of this paper we focus on the leading order 
free energy foits) — limAr^oo F{ts)- Expectation values of bubble observables have similar expansions. If B is 
a D-colored graph, the connected expectation value 



1 \TrB(T,T); _ 1 ^^.(^) dF^ y (-1)H Upte,,, ^2.13) 



N Z ND dtB ^ s(g) 
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Figure 3. Trace invariants and gluings of simplices in D = 3. 



has an expansion in connected {D + l)-colorcd vacuum graphs Q having S as a (marked) subgraph, denoted G Z) B. 
The scaUng in N in (|2.13p of a graph Q rewrites 

Ar-j5T"(e)7v~^5u^("(^^"-°"'^') . (2.14) 

Using (a weaker version of) proposition [2] in appendix!^ ^{G) ^ Dijj{B) and the inequahty is saturated for io{G) = 0. 
It follows that in the large N limit only graphs Q D B of degree zero contribute to the expectation. 



C. Topology from bubbles 



To simplify the discussion, in this section we will restrict to the case = 3. The original idea of tensor models 
[28l - l30j was to generate triangulations of 3-dimensional spaces. The basic building block in the original proposals was 
an interaction term which combinatorially describes a tetrahedron (a 3-simplex) also used in group field theories (ssj 

Vtotiahodron = TabcTcdeTebfTfda ■ (2-15) 

This term is not (E)^U{N) invariant. The most one can say about it is that it is invariant under a simultaneous 0{N) 
orthogonal transformation of all its indices. 

The situation is already improved in colored tensor models [s^l where the indices are distinguished and one can 
implement a 0^U{N) invariance. As the pattern of contraction of a tetrahedron is not a trace invariant one can raise 
the question of the topological interpretation of the trace invariant observables and their relation to triangulations. 

The situation is actually like in one-matrix models with generic interactions. A Tr(Af '^)-vertex is seen (by duality) 
as a polygon with k sides. A closed graph is then a gluing of such polygons. Obviously one can divide each polygon 
into triangles (by adding a vertex in the middle of the polygon, i.e. by taking the topological cone over its boundary), 
so that the graph encodes a triangulation. Here, a similar interpretation holds. The (3 + l)-colored graphs are known 
to describe topological 3-dimensional pseudo- manifolds [s^l- The black and white vertices of the graph correspond 
to tetrahedra (3-simplices). The triangles (2-simplices) bounding a tetrahedron are represented by the half-lines 
touching the vertex, hence are colored 0,1,2,3. The lower dimensional simplices are colored by the colors of the 
triangles sharing them. Thus the edges are labeled by pairs of colors (the edge 12 is common to the triangles 1 and 
2) , and the points (vertices of the tetrahedra, to be distinguished from the vertices of the graph) are labeled by triples 
of colors (the point 123 is the point common to the triangles 1, 2 and 3 bounding a tetrahedron). 

A line in the colored graph represents the unique gluing of two tetrahedra of opposite orientations along boundary 
triangles which respects all the colorings: that is we glue triangles of the same color, say 2, in such a way that the 
edge 02 (resp. 12 and 32) bounding a triangle is glued on the edge 02 (resp. 12 and 32) bounding the second triangle, 
and similarly for points. This construction yields the pseudo-manifold dual to a (3 -I- l)-colored graph. 

Alternatively the same graph with 3 + 1 colors can be seen as the gluing of the effective interactions, B which are 
graphs with 3 colors, along the lines of color 0. Following the above construction, each effective interaction by itself, 
being a graph with 3 colors, represents a surface. The (black and white) vertices are dual to triangles, and the edges 
bounding the triangles are colored 1, 2 and 3. The surface represented by an interaction is the unique one obtained 
by gluing the triangles along their edges (as indicated by the graph with three colors) respecting all the colorings (i.e. 



those of the edges and of the points). In figure 3(a) for instance we represented such a surface obtained by gluing 
eight triangles. 

Adding the lines of color results in taking the topological cone over this pseudo- manifold. Cm = {M x [0,1])/ (M x 
{1}). Let us first examine the effect of this coning on one triangle (represented in figure [3(b)| ). The original triangle 



(a) (b) 
Figure 4. Stranded graph and a non-bipartite invariant. 

will now be called a triangle of color (sec figure [3(b)| ) . The original edges acquire the new color 0, hence they will be 
called 01, 02 and 03 (see again [3(b)| ) , and similarly the original points (13 becomes 013, etc.). This coning adds extra 
triangles, edges and points. Every original edge gives by coning a new triangle. We color this triangle by the color 
of the edge, hence the original edge of color 1 gives rise by coning to the triangle of color 1 (see again figure |3(b)| . 
Note that the new triangle 1 shares with the original triangle, the edge 01. Every original point gives by coning an 
edge, which inherits the colors of the original point (the edge 13 is the cone over the original point 13 and is shared 
by the triangles 1 and 3). We also obtain a new point labeled 123. When taking the cone over the surface defined by 
a connected graph B (with colors 1, 2 and 3), we obtain new triangles (one for every edge of the surface), new edges 
(one for every point of the surface), and an unique new point 123, the apex of the cone. 

Thus, when seen as a subgraph of a (3 + l)-colored graph Q, B represents a "chunk" of the 3-dimensional space. 



For the example of the graph with 3 colors in figure 3(a) adding the dashed lines of color 0, we obtain the gluing of 



8 tetrahedra drawn in figure 3(c) This chunk has the topology of a ball and is bounded by the 8 triangles of color 
corresponding to the dashed half-lines. 

Thus the trace invariant quartic interactions like 



represented in figure 
0, and not to a tetra. 
torus. 



3(d)[ correspond to a gluing of four tetrahedra, with four external, boundary triangles of color 



redron. Note that a chunk can have a non-trivial topology, for instance it can be a cone over a 
One can employ an alternative stranded graph representation of the Feynman graphs, closer to the ribbon graph 



representation of matrix models. This is presented in figure 4(a) One replaces the black and white vertices of the 
effective interactions by stranded half-lines, which are then connected by dashed lines having each three strands. In 
this representation the strands colored 1, 2 and 3 have each an associated Kronecker 5 which corresponds to the 
contraction of a tensor index between two tensors of the bubble observable. The dashed lines have three strands 
representing the three Kronecker 5 coming from a Wick contraction which propagate the tensor indices. The faces 
of colors Oi arc easily identifiable. Each stranded half-line corresponds to a triangle (the triangles of color in the 
colored graph representation). The graph of the effective interaction encodes the pattern of gluing of the triangles 
into a surface (the boundary of a chunk), and the dashed lines encode the gluing of chunks along boundary triangles. 
As this representation is redundant and somewhat cumbersome we will not use it further. 

Before concluding this section let us remark that the fact that the graphs are bipartite plays a secondary role, 
ensuring just the orientability [6^ . What is crucial is that a colored graph represents the unique gluing of simpliccs 
which respects all the labellings (including the induced ones over all the lower dimensional simplices). Dropping 
the bipartite requirement allows one to consider the 0{N)^ invariant presented in figure |4(b)[ As it consists in a 
gluing of four triangles and any two triangles share exactly one edge one might be tempted to interpret it as a gluing 
pattern of four triangles bounding a tetrahedron. However, this interpretation is not the correct one. Indeed, on 
closer inspection, it turns out that the dual gluing consists in four triangles glued first around a vertex (say 13) and 
then glued along opposite edges of color 2 (sec the right hand side of figure [4(b)| ) . Thus, respecting the rules of the 
colored gluings described above, this 0{N)^ invariant has the topology of the real projective plane RP^. 

III. LARGE N LIMIT 
A. The melonic family and the large N factorization 



In the large N limit, only graphs with vanishing degree survive. For (2 -I- l)-colored graphs the degree is the 
genus of the graph, hence the graphs of degree are exactly the planar graphs and represent spheres. For I? > 3, 
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the (D 

Mm 



l)-colored graphs Q with uj{G) = have been shown to also describe topologieal spheres in dimension D 



Combinatorial description of melons 



We explain in appendix lAl why the (D + l)-colored graphs of degree zero, called melonic, can be obtained by the 
insertion procedure detailed below. While the dominant graphs of our models are melonic, it is understood that not 
all melonic graphs are generated. However, this section is only concerned with the combinatorial properties of the 
melonic family, hence we temporarily allow ourselves to also use melonic graphs which do not appear in the Feynman 
expansion of our models. 

First order. The lowest order graph consists in two vertices connected hy D + 1 lines, as in the figure 5(a) We 



consider all lines incident at the positive (black) vertex to be active, which means that higher order graphs will be 
obtained by insertions on them. 

Second order. D + 1 graphs contribute to the second order. They arise from inserting two vertices connected 
by D lines on any of the D + 1 active lines of the first order graph. If the line on which we insert this decoration 
has color i, the new lines will be colored by all colors except i. Say we insert this graph on the active line of color 1, 
like in the figure [5(b)] (hence the new lines have colors 2, 3 up to D). All lines incident at the new black vertex are 
deemed active (the new lines of colors 2, . . . ,D as well as the external line of color 1), while the exterior line of color 



1 incident at the new white vertex (in bold in figure 5(b) ) is deemed inactive. 




(a)Tho first order melonic 
graph. 



(b)A second order melonic 
graph. 



Figure 5. 

Order p+1. We obtain the graphs at order p + 1 by inserting two vertices connected by D lines (with appropriate 
colors) on any of the active lines of a graph at order p. Once again, with respect to the new vertices, all lines incident 
to the black vertex are deemed active, while the exterior line incident to its white vertex is deemed inactive. 

We are now going to show that the expectation values of melonic graphs are fully determined by the (dressed) 
covariance of the model, in a specific, factorized form. 



2. Large N factorization 



In the large N limit, only the bubble observables B for which there exist (D + l)-colored graphs Q which are melonic 
survive. The melonic graphs have some important properties, which put together lead to the large N factorization of 
expectations. 

• If a (D + l)-colored graph Q is melonic then all its subgraphs B with colors 1,2, ... ,D are melonic see figure 



6(a) and are therefore built following the same procedure. 



• In this procedure, Q is obtained by inserting pairs of vertices v and v separated by D lines, and a D-colored 
subgraph B is obtained by performing the same insertions, but ignoring the color 0. 

• Consider two vertices v and v inserted at some step. At the time of the insertion, they are connected in Q by 
D lines and some two-point graph (corresponding to the line on which they have been inserted). As all further 
insertions arc made on the lines of Q, the two half-lines of any color (0, 1, up to D) on v and v will always be 
connected together via two-point graphs. 

Therefore, for every such pair of vertices v and v of B, the two half-lines of color must be connected via some 
two-point graph in Q (see figure [6(b) ). In other words, starting with a melonic D-colored observable, there is a unique 
way to pair its external half-lines with two-point insertions so as to get melonic (D + l)-colored graphs. Then, the 
full expectation value is obtained by inserting this way full two-point functions, one for each pair of vertices which 
are joined. 



2 

(a)Melonic observables. 



(b) Graphs contributing to the melonic 
observables. 



^1 



Figure 6. Graphs contributing to the 3-dipole expectation. 



As a result , in the — > 00 limit , the expectation of a melonic observable factors in terms of full two-point functions 
(dressed propagators). The full two-point function writes 



T T- - 



(3.17) 



where j- is the bare propagator and the dots denote the radiative corrections. We denote limAr^oo U (tg, N) ~ U [te). 
The large expectation of the ZJ-dipolc observable Bi computes then 



lim 



1 (Tr8,(r,T)) 



= U{tB) 



(3.18) 



N^oo N Z 

Some graphs contributing to this expectation for 13 = 3 are presented in figure [6l where the marked graph Bi is 
presented in bold. The colors of the lines are assigned turning clockwise 0, 1, 2 and 3 at the white vertices. 
We thus conclude that 



1 \TrB(r,T) 



Z 




if B is not melonic 

if B is melonic with 2pg vertices 



(3.19) 



In particular this factorization holds for the Gaussian model with = for B ^ Bi, and the full two-point function 
simply given by the bare covariance [/(tg) = 1/ti. 

The universality of tensor measures, first derived in [56j (where it was obtained by mapping melons to trees), is the 
fact that the observables satisfy (|3.19p . It means that in the large A^ limit the models become Gaussian. However this 
large A^ limit is very non-trivial, as the covariance of the large A^ Gaussian is the full, resummed, two-point function. 
The rest of this paper is dedicated to explore the various multicritical behaviors and continuum limits governed by 
this resummed covariance. 



B. The leading order two-point function and free energy 



The full two-point function at large A^ is determined by a self-consistency equation provided by a Schwinger-Dyson 
equation supplemented with the above factorization. The relevant Schwinger-Dyson equation is 



— y - 

ni....,no 



dTdT 



d 



rp -N"-^S{T,T} 

J- TLi. ..nr) ^ 



. 



(3.20) 



Taking the derivative explicitly, one gets 



Tr6(r,r) 



= 0, 



(3.21) 
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where ps denotes the half- number of vertices of the bubble B (that is the number of either black or white vertices). 
At leading order in 1 /N this can be rewritten in the following form. We first define the leading order potential 

F(x,%) = E( E ^e)^"' ^'(x,tB)^|^(x,tf5) = 5]n( (3.22) 

n>l B mclonic, n>l Z3 mclonic, 

PB=n PB=n 

and taking into account the factorization of the melonic expectations, the Schwinger-Dyson equation becomes the 
following self-consistency equation 

UitB)V'{U{tB),tB)^l . (3.23) 

The leading order two-point function is the solution of this polynomial equation whose coefficients are the coupling 
constants of mclonic obscrvablcs. 

Once U is determined using (|3.23p . one can access the free energy /q. The leading order free energy foitg), like 

the leading order potential V{x,ti3) — X]n>i(X]s mclonic. ^s)^;" and the leading order two-point function C/(te) only 

depends on the coupling constants of the melonic bubbles tjs. Consider the function /o — V[U {tB)TtB) + \nU (tg). Its 
differential is 



d\fo~V{UitB),tB)+lnUits)]= \ji^-U{tBr ^V\U{tB),tB)^ + ——^]dts^O. (3.24) 

o mclonic 

Thus the leading order free energy is 

k{tB) = V{U{te),t,s) - In [/(te) . (3.25) 



C. The continuum limit 



Tensor models arc of combinatorial nature and as such provide a notion of continuum limit in a combinatorial way. 
The idea is that disregarding the geometrical content and interpretation which may be given to a model, this limit 
is always obtained as the regime where graphs with a very large number of vertices dominate. As an illustration, we 
derive the continuum limit in a particular truncation defined by the action 

5*7-4 (T, f) = ^ T,tf,-J + 5 ^ Taj^...ao-iaDTai---ao-ibD'^bi---bD-ibn'^bi---bo-iaD ■ (3.26) 
fl ai,...,ar) 
bi,...,bo 

Note that the interaction term is melonic. The leading order potential is defined by ti = 1,^2 = .9, that is V{x) = 
X + gx^ . The leading order free energy is therefore 

/o(5) = E ff"/'""' ' (3.27) 

nelN 

where /(^") is the number of {D + l)-colored melonic graphs built with n effective interactions (thus having 4n 
black and white vertices). The number /(^") is a canonical partition function for graphs with fixed number of vertices 
and fo{g) is its associated grand-canonical partition function with lattice "chemical potential" g. The thermodynamic 
limit is encoded into the asymptotic behavior of /'■''"•', 

/(4») _ A n''-'^ 5"" , (3.28) 

for some constants g^ and 7. Thus gc is the radius of convergence of foig), which means that when g approaches 
9c, foig) loses its summability and graphs with a large number of vertices (4n) dominate its behavior. The power-law 
decay characterized by 7 controls the singularity of fo{g) close to gc, since 

Mg) - \g ~ gd'--' . (3.29) 

The exponent 7 is known as the entropy exponenll^. Let us compute the entropy exponent of the model defined 
by S't*. First one notices that the derivative of the leading order free energy writes in terms of the leading order 



^ In the Tr Af* matrix model for random two-dimensional lattices for instance one has 7 = -1/2 Q, which is the universality class of 
pure 2d quantum gravity. As such, it is reached generically, i.e. for most values of the coupling constants, in the continuum limit of 
one-matrix models. 
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two-point function = U{g)'^. The equation (|3.23|) gives 



2gU{gf + U{g)- 1 = 0, hence Uig) = , (3.30) 

where we selected the physical root with initial condition C/(0) = 1. One thus identifies gc = —1/8 and for g gc the 
non-analytic parts of the two-point function and free energy are 

U (g)sing - (g -gc)K /o,si„g - (g - gcf^^ , hence 7 = 1/2. (3.31) 
The average number of effective interactions (proportional to the number of vertices) diverges when tuning to criticality 

d 1 

{n) = g ^ log /o - : 00 . (3.32) 

dg \g-gc\ 9^9a 



D. Virasoro constraints 



We now come back to the generic model with arbitrary couplings. The full set of Schwinger-Dyson equations is 
obtained by inserting generic _D-bubble observables in (|3.20p . The equations can be recast as ig2'(tg) = for some 
differential operators Lb labeled by the observables. The algebra of these operators has been discussed at length in 
[soj . Due to the large N factorization one can find the leading order Schwinger Dyson equations and the associated 
algebra of constraints by a shorter route. We show below that the large N factorization (|3.19p reduces the Schwinger- 
Dyson equations to a set of Virasoro constraints, like in matrix models. We emphasize that this only holds at leading 
order in 1/N . 

Note that in fact the leading order two-point function U{tB) as well as the leading order free energy /o depend 
only on the sums of the coupling constants of melonic observables at fixed number of vertices. Thus, defining 

t„ = melonic. ^B, thc large N factorization becomes 

Pe=n 

^ = Pit,)]'' , (3.33) 
Then multiplying (|3.23p by U'' for any positive k we get 

[Uitp)f - ^ n i„ [C/(tp)]"+^ = , Vfc > . (3.34) 

ri>l 

These equations can be seen as differential equations on /o, 

ifc /o = , for Lk = tt: " tn tt: ■ (3.35) 

dtk ^ Otn+k 
n>l 

It is straightforward to check that the differential operators satisfy the well-known algebra 

[L„, Lrn] = {m-n) Lm+n ■ (3.36) 



IV. MULTICRITICAL BEHAVIORS IN THE LARGE N LIMIT 



We now set the summed coupling constants ti = l/g, and tn = an/g to investigate the different possible critical 
behaviors with respect to the parameter g. The potential V becomes V{x) = {x + J2k=2 x^)lg- The self-consistency 
equation (|3.23p for the two-point function becomes 

m 

g = C/ + ^fcafc[/^ (4.37) 

k=2 

Like in the matrix models review [l|, the continuum limit is obtained by setting dg/dU = 0. At least locally, this 
equation can be solved for U in terms of the parameters a^. One concludes that for generic values of these parameters 
7 = 1/2, thus proving the universality of this continuum limit, first derived in (4l| for the colored model. 
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But there are points on the set of parameters where the equation dg/dU = can not be solved for U, because 
d^g/dU^ may vanish and there the impHcit function theorem can not be apphed. In these cases, one has multicritical 
behaviors and 7 > 1/2. A multicritical point of order m is defined, like in one-matrix models, by 



dg_ 

dU 



d 



im— 1 , 



0, 



and 



9™5 



^0 



(4.38) 



which imply g = gc — {Uc — U)"^ + 0((U — C/c)™"''^). Such multicritical behaviors have already been observed in 
tensor models in [48| (where they are interpreted in terms of dimcr models) and 53|. As 17 is a polynomial in U whose 
coefficients a^. can be freely chosen, multicritical points can be reached for the generic one-tensor model. We present 
below a minimal realization, i.e. a potential V with minimal degree leading to a multicritical point of order to. We 
first set Uc = to^ , and consider 



1 . . 1 

3 k=l ^ 

Thus V is of degree to, satisfies V{Q) ~ 0, the coefficient of the linear term is 1/ g and UV'{U) — [U^ 
The self-consistency equation is exactly 



g = gc-{Uc- U) 

Substituting into p.25p we find 
1 ™ 1 



with 



9c = u: 



(4.39) 

{Uc~Uy"]/g. 
(4.40) 



U'-iUr-U)" 



fo{gc) 



q '■ — ' k 

^ fe=l 



-\^{Uc-{Uc-U)\ 

k 



E 

fe=i 



f fUc-U 



k 



Uc 



1 1 

9 9c 



m ^ 



(441) 



k=l 



Taking into account that Uc ^ U ~ (gc — g) , we obtain for g — > gc 



/o = /o(.gc)+(Ei 



k=l 



9c- 9 



(to + 1) 



1+^ 



+ Oiigc-g)''-^) 



(4.42) 



hence a multicritical entropy ex pon ent 7™ = ^ — a-s obtained in [5^ and in [l^, coinciding with the ones of 
multicritical branched polymers [55| . 
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Appendix A: Combinatorics of colored graphs: jackets, degree and melons 

In order to define the degree of a graph one first needs the notion of jacket (ssl-lioj . 

Definition 2. Let B be a D-colored graph and t be a cycle on {1, . . . ,-D}. A coZorec? jacket J of B is a ribbon graph 
having all the vertices and all the lines of B, but only the faces with colors (r'(l), r''^^ (1)), for q = 0,...,-D — 1, 
modulo the orientation of the cycle. 

As a jacket J oi B contains all the vertices and all the lines of B, J and B have the same connectivity. As such, any 
jacket J carries some key topological information about B (for instance the fundamental group of S is a subgroup 
of the fundamental group of any of its jackets [13]). For graphs with four colors, the jackets correspond to Heegaard 
splitting surfaces [49|. 

Jackets are ribbon graphs, hence they arc completely classified by their genus gj. 

Definition 3. The degree Lij{B>) of a colored graph B is the sum of genera of its jackets, id{B) — gj- 
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Graphs with three colors are ribbon graphs, and the degree coincides with the genus. The crucial property of the 
degree is that the total number of faces of a D-colored graph computes in terms of the degree. 

Proposition 1. Let B a D-colored graph with 2p vertices. Then the total number of faces of B respects 

Proof: This equation can be found in the literature (see for instance). However, due to its importance we 
present here its proof. Every jacket is a ribbon graph with 2p vertices and Dp lines, hence the number of faces of a 
jacket is 

\J^j\ = {D^2)p + 2-2gj . (A.2) 

As jackets correspond to cycles over D elements modulo the orientation, B has ■^{D — 1)! distinct jackets. The faces 
with colors ij will belong to the 2{D — 2)1 jackets corresponding to the cycles tt such that 7r(i) = j and 7r(j) = i. 
Moding by the orientation we conclude that each face belongs to exactly (Z) — 2)! distinct jackets. Summing (jA.2[) 
over the jackets and dividing by ^{D — 1)! proves the lemma. 

□ 

Of course the same definition goes trough for graphs Q with D + I colors. Further facts concerning the degree are 
listed below. 

Proposition 2. The degrees of a {D + l)-colored graph Q and of its D-bubbles B^p) with colors 1, . . . ,_D respect 

L^iG) > Dj2^iB^p^) . {A.3) 

p 

The proof of this statement can be found in [l^], lemma 7. The proof relies on the identification of the jackets of 
S(p) as ribbon subgraphs of the jackets of G- 

Proposition 3. // the degree of a [D + \)-colored graph Q vanishes, i.e. all its jackets are planar, then Q is dual to 
a D-sphere. 

The proof of this lemma can be found in [i^ . 



The graphs of degree have been thoroughly analyzed in [41[. Their characterization relies on two lemmas. 

Proposition 4. If D > 3 and Q is a (D + \)-colored graph with vanishing degree, then Q has a face with exactly two 
vertices. 

Proof: All faces of Q have an even number of vertices. Denote \Fs\ the number of faces with 2s vertices. By 
proposition [1] (taking into account that Q has D + 1 colors), the total number of faces of Q is = X]s>i = 

^^^2~^^ P + As a vertex belongs to ^'^^^+^'> faces we have X]s>i ^l-^s| ^ ^^^2^^^ P- Eliminating J^2 we get 

:F, = 2D + Y^is 2)Fs + W-^ p . (A.4) 

s>3 

The first two terms give a strictly positive contribution for any D, whereas the third term changes sign for D ~ 3. 
Thus J"i > 1 only for L> > 3. 

□ 

The fact that this proposition fails in D = 2 is the source of the difference between the large limit of matrix 
models, dominated by planar graphs, and the large N limit of tensors of rank D, dominated by melonic graphs, which 
we describe below. 

Proposition 5. If D > S and Q is a {D + l)-colored graph of vanishing degree, then it contains two vertices v and v 
separated by exactly D lines. 

The proof of this lemma can be found in j4l| . It relies on proposition |4l but as it is somewhat convoluted we do 
not reproduce it here. 

We exploit this lemma in the following way. Starting from a {D + l)-colored graph Q we identify two vertices v and 
v separated by D lines. Erasing this subgraph and reconnecting its external lines we obtain a graph having two less 
vertices and degree (as it can be checked explicitly). Iterating this erasing procedure we necessarily end up with 
a (Z? + l)-dipole, that is the graph having two vertices connected by I? + 1 lines. Conversely, every melonic graph 
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can then be obtained by starting from the {D + l)-dipole and inserting such subgraphs (consisting in two vertices 
connected by D hues) arbitrarily on all the lines. 
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